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ABSTRACT 


Recent  suggestions  have  been  made  for  the  design  of  bulk  -semiconductor 
millimeter  wave  devices,  and  in  particular,  a  new  type  of  phase  shifter. 

In  order  for  these  designs  to  perform  in  a  satisfactory  manner  it  was  found 
necessary  so  demonstrate  that  electromagnetic  propagation  would  occur 
largely  in  the  interior  of  a  semiconductor  dielectric  waveguide  with 
relatively  low  loss.  In  this  report  an  analysis  is  made  of  electromagnetic 
waves  guided  in  an  infinite  slab  of  material  with  the  properties  of  high 
resistivity  silicon.  Calculations  end  preliminary  experiments  are  demonstrated 
at  frequencies  near  16,0  GHz.  It  is  concluded  that  propagation  in  the  semi¬ 
conductor  medium  offers  the  possibility  of  low  loss  circuitry  and  satisfies 
the  requisites  for  device  design.  , 
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BULK  SUII CONDUCTOR  QUASI- OPTICAL  CONCEPT  FOR  GUIDED  WAVES 
FOR  ADVANCED  MILLIMETER  WAVE  DEVICES 


INTRODUCTION 

In  a  recent  proposal^-  we  have  suggested  the  design  of  bulk  semiconductor 
devices  end  the  use  of  semiconductor  waveguides  to  replace  more  conventional 
stripline  and  metal  waveguides  for  the  millimeter  wave  regions  of  the 
spectrum. 

It  is  well  known  that  metallic  waveguides  become  lossy  at  shorter 
wavelengths  (1  dB  per  foot  at  3*2  mm  or  1.5  dB  per  foot  at  2,15  mm).  The  use 
of  dielectric  guides  offers  a  promise  of  lower  losses  as  the  wavelengths 
approach  the  optical  regions.  With  the  advent  of  recently  obtainable  high 
resistivity  silicon,  low  loss  semiconductors  offer  an  attractive  potential  89 
a  guide  material.  The  additional  factor  here  is  that  components  and  devices 
could  conceivably  be  "built  in"  the  semiconductor  line  so  that  in  effect  one 
would  have  an  optical  line  with  active  elements  incorporated.  The  first 
considerations  Involved  in  evaluating  such  a  concept  is  the  direct  test  of  a 
single  straight  section  of  semiconductor  dielectric  guide  to  determine  the 
losses  which  might  result.  A  factor  in  this  line  of  thought  is  that  if  the 
energy  flow  is  almost  entirely  within  the  waveguide,  it  would  be  relatively 
easy  to  control  by  electrical  means.  If  the  energy  flow  is  along  the  direc¬ 
tion  of  the  guide  but  contained  in  the  air  space  surrounding  the  semiconduc¬ 
tor  guide  it  would  be  relatively  difficult  to  control.  See  Figure  1, 

The  firBt  consideration  is  to  estimate  the  chances  of  power  flow  within 
the  guide  material,  where  it  can  be  modulated,  amplified,  detected,  etc.,  by 
devices  constructed  and  located  within  the  semiconductor  structure.  In  what 
follows  we  Bhall  describe  first  some  analytical  considerations  of  the  semi¬ 
conductor  system  relying  on  the  wave  guided  in  a  dielectric  slab.  This  will 
then  be  followed  by  some  initial  experiments  designed  to  prove  that  almost  all 
of  the  energy  is  contained  in  the  semiconductor. 

DIELECTRIC  SIAB  WAVEGUIDE 

We  consider  a  dielectric  slab  oriented  in  the  Indicated  position  in 
Figure  2.  Here  we  assume  an  infinite  dielectric  slab  with  no  variations  of 
field  in  the  x  direction. 


O 

The  general  form  of  Maxwell's  Equations  are  : 
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i l£  ±  H  j.T£  -  -  1 

'  '  i>  (a) 

£j± +  ££  +  ££.-.  v^eii  J 

**'  >y3  TP'  ' 


Propagation  is  in  the  z  direction  so  that  we  use  the  form 


<r  +  ^/3. 


Purthermore  E  end  H  can  be  a  function  of  y  but  cannot  be  a  function  of 
x,  bu  assumption.  This  means  that  all  terms  with-E—  =  0,  and-fl-  is  re¬ 
placed  by  —  Y,  »*■ 


We  next  take  this  into  account  in  (1). 
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These  equations  can  now  lead  us  into  the  calculations  of  either  TM  or 
TE  modes.  We  choose  the  former  to  illustrate  physical  principles. 
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We  also  require 


-  -  u//t  £l 

n;=-^n; 


Next  we  consider  the  fields  outside  the  dielectric  slab,  y} 
and  inside  the  dielectric  slab  y< 

First  consider  the  fields  outside  the  slab,  y  } 


V *H,  -  ~  <r  /f), 

&  It*  ~f-  Y&  /fy  ~  (r  //y 
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Note  we  have  the  following  terms  as  a  function  of  Hx 


-  c  -? '  k’y 

£y  Hj 

r-°  e 

CO  6  tr 


Lx  - 


*  Si 
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df  P 
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This  is  related  to  rectangular  guide  as  in  Figure  4. 

Inside  the  slnb;^/^  Y,  }  % 

In  the  z  direction  of  propagation 

Hi  -  H^er,i  In) 


Hoter^r  Y0  in  order  for  fields  to  match  across  the  surface. 

Also  by 
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b* 
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-  k)  a/ 
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»/ 
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These  tvo  cases  are  indicated  as  in  Figure  5. 


The  surface  impedance  is 


-  -Jz*.  -  jJx  . 
5  /4  U(l 


Since  hQ  is  realf  ^  is  inductive  reactive. 
In  summary  we  have  the  equations  shown  In  Table  I. 


Table  I:  Equations  for  the  TM  Mode 


Outside  the  Slab 

i>/r  h?Hx° 


Inside  the  Slab 


/V/-  Ce'A,/ 
-  -jjiy* 


H,  - 

=  -  k,1  h/ 

tr  kf  r  JjJ  'b  6, 

~  stn  h,  y  Ant  CoS  h,  y 

V  -  * 

’  £g  =  ^  XMf 

=  jJh  c'a>  h,y 
u>  6,  dy  u>  6,  ' 

or  £  -r^_  h,y  -  -J_hj  5-//,  4  !/ 

t,  -  * 
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Consider  two  possible  modes  of  propagation. 

Case  I  Jbr  -a/2  (  y^a/2,  assuming  =  sin  k^v 

H /.  ~  B Q  •>  >r\  h, y ^ 

H,-Ack\  y>%  ■ 

At  the  boundary,  y  =  a/2 

Ac  8s, *(*>,%).  (jx j 

Also  at  the  boundary 

~  by  (3*). 

In  general, 

/ 2  - ^ -  Be  *s,«A,y 

^  by  A 

E*.  =  -&Bc'~'hY  •  -#</<*> 


(?■// 


for  y  ^  a/2 
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Hence  at  y  =  a/2  f 

-t-% 


/?e 

-V  ,  i  ^ 


—  jii — -  8  k,  ^ ,  (zj J 

fu>t,  -A 


By(’l8jand(23^ 
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±j>< 


i — la.  — 


h. 


-  - 
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M 


In  addition: 


C— 


v  = 


t- 

r+ 


tr 

“!>'  i, 


by  definition 


adding 


/  ^  /  A  t 
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<?</ 
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Multiplying  (24)  by  a/2,  and  (26)  by 

c4 

-  ^  £  Cert  4i)  - 

4  tL. 

47) 

A  V 

-hi  cl.  —  10 s*f(,  ' 

■  <>) 

r^j 

> 

> 

II 

/  = 

(?0 

c,+f  - 

C 3c J 

) -!  : 

-  /r2- 

(7) 

y 

For  Ku  Band, 

a  =  .3U  x  2.54  =  7.9xl0_1cm  =  7-9  x  10”3!!! 

u)  «  2 trf  =  2  7T  x  15  x  105  =  94  x  109  =  9.4xl010 

6,  =  12  x  8.854  x  10"12  =  106  x  lO'22  =  1.06  x  10"10 

£  =  1  x  8.854  x  10"12  =  8.854  x  10"12 

ir 

L  =  1.257  X  10" 6 

0,  =  11  x  8.854  X  10’12  =  97  X  10-12  =  9.7  X  10mU- 
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^  -  (9.4)Scl020  x  1.257x10’ 6  x  9.7xlO"n  x  (7.9)2  x  lo’6 

4  - 

10  o 

R  =  9.4x10  x  fTT257  x  10  x  x  lo”6  x  7. 9x10" 3 

- 2 - 

=  9.4  x  ^257  x  ffi  x  7.9  x  1010  x  lo’3  x  10*6  x  lo’3 

~T 

3  9.4  x  1.12  x  9.9  X  3.95  X  lo'2 


R  =  4.4 
R2  =  19.3  . 

A  complete  plot  of  30  and  31  vas  inm  solving  for  p  and  q. 
The  results  are  Indicated  In  Figure  6. 

•From  this 
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/ 

/y/ 


-3-  , 

P  7^b n 


The  roots  for  equations  (1»0)  and  (1*3)  are  shovn  in  Figure  8. 

Case  II  a. 

P  ~  /.6'SfC 

e  -  y.  /«*■  9-  Wi) 


Solve  for 


*&/£*■ /.&  -  Ainu/ 

7fX,°  3 

,  y  =  /°'s 


c.  e. 


/  Shf  X 3-  -  ,  £Cj  X/oJ 
7  ?  too-' 

CoS ' J.<  ~  ccs  /  S  3 ft  ~  .  Cioi 
5+e 
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Case  II  b. 


A 

V'3/5  y  %  ,  /d  S 

A-  -  S4- 

—  /■  7-  xjc  *  —  .2.  aj  y /o 

■V =/'  , 

y  /o'*' 

J./y 

(^'/J 

k, 

T-  X,  t  j_  Zjd  ^  /  09  X  /O^ 

?■  9 

(**) 

C°s(kt 

—  —  Co**/- 1'*'—  — 

This  tells  us  that  in  this  mode,  more  of  the  field  is  external  to 
electric  and  hence  there  is  a  greater  chance  for  loss. 

the  di- 

EXPERIMENTAL  DATA 

A  laboratory  bench  was  set  up  with  the  components  arranged  as  indicated 

in  Figure  12.  This  is  essentially  a  bridge  arrangement  with  a  klystron 
source  capable  of  supplying  Ku-band  radiation,  a  bridge  with  an  attenuator 
and  phase  shifter  in  one  arm  and  sections  of  waveguide  which  can  be  removed 
in  the  other  arm.  A  standard  diode  detector  was  used,  the  output  of  which 
was  fed  into  an  oscilloscope.  The  tuner  is  also  indicated  in  the  figure. 
The  section  of  guide  labelled  "A"  could  be  removed  and  four  post  conditions 
were  studied  as  indicated  in  Figure  13.  In  Case  A,  the  empty  waveguide  is 
filled  with  silicon  ^7/32"  long  with  cross  sectional,  area  of  .311  x  622", 
l.e. ,  fitted  into  the  inner  dimensions  of  a  Ku-band  waveguide. 

The  following  procedures  were  used  particularly  in  the  cases  of  Condi¬ 
tions  B  and  C. 


1.  The  attenuator  was  set  at  maximum,  i.e.,  greater  than  50  dB 
so  that  all  of  the  power  would  flow  through  the  am  containing  the  semicon¬ 
ductor. 


2.  The  klystron  frequency  was  changed  slightly,  at  the  same  time 
adjusting  the  tuner  on  the  detector  for  maximum  power  transfer.  Usually 
three  or  four  peaks  were  found  which  were  indications  that  at  these  perticu- 
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lar  frequencies  the  semidonductor  block  was  acting  as  a  l/2  wavelength 
transformer.  It  was  found  that  these  peakB  were  approximately  .78  GHz  apart. 
The  data  obtained  1 «  nViovn  in  Table  II  indicating  the  EF  frequency  at  test, 
attenuation,  and  changes  in  attenuation  over  the  control.  From  this  table, 
it  can  be  seen  that  the  silicon  filled  waveguide  gave  an  added  attenuation  of 

4  dE,  while  the  silicon  waveguide  with  no  cover  gave  an  added  attenuation  of 

5  dB.  We  conclude  from  this  that  one  additional  dB  of  attenuation  was  found 
by  completely  removing  the  metal  cover  from  the  silicon. 

TABLE  II 


TRANSMISSION  DATA 


Condition 

frequency 

Attenuation 

Attenuation 

gh2 

dB 

AControl  -  ATest 
dB 

Air  filled 

(A)  waveguide 

(control) 

16.50 

9-7 

9.2 

9.0 

9.5 

5.35  Average 

0 

Silicon  filled 

15.40 

13.0 

(B)  waveguide 

15-46 

13.0 

4 

16.18 

13.5 

Silicon  waveguide  15. 54  14.1 


13.9 

2 5-"  exposed  l6,40  lU.y  5 

14.5 

(C)  to  air  17.16  14,6 

Silicon  removed 

2i"  air  space  17.04  26.5  17 

0» 


We  should  like  to  discuss  the  interpretation  of  the  data  in  Table  II. 

1.  The  control  test  showed  no  resonance  points  under  Condition  A. 
There  was  no  sign  of  l/2  wavelength  transformer  action  as  might  be  expected. 

2.  Under  Condition  B,  l/2  wavelength  transformer  action  was  ob¬ 
served.  The  peaks  in  transmission  are  assumed  to  eliminate  surface  reflec¬ 
tions  and  the  4  dB  attentuation  which  was  noted  can  be  attributed  entirely  to 
the  losses  due  to  conductivity  in  the  silicon. 
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3.  Under  Condition  C,  measurements  were  made  at  peak  transmission 
values  for  the  same  reason  the  attenuation  was  found  to  he  5  8B  indicating 
additional  1  dB  loss.  The  exposed  silicon  vaa  probed  with  metal  stmeturen 
such  as  a  screvdriver  or  metal  plate.  Little  or  no  effect  on  transmission 
was  found  by  placing  flat  plates  on  the  top  of  the  silicon  guide  or  on  the 
sides  of  the  silicon  guide  near  its  center  of  exposure.  However,  when 
placing  a  round  metal  rod  approximately  l/8"  diameter  near  the  corner  of  the 
exit  of  the  waveguide  flange  and  the  silicon  on  the  upper  flat  piece,  a 
slight  change  in  balance  was  noted  of  approximately  0,2  to  0.3  dB.  Changes 
from  this  preliminary  probing  appear  to  show  that  there  may  be  some  indica¬ 
tion  of  power  loss  due  to  a  poor  fit  between  the  silicon  and  the  waveguide 
wall.  This  must  be  further  examined  and  if  bo,  there  are  still  two  approaches 
to  overcome  this  problem.  First  a  horn  can  be  used  to  launch  the  wave. 

Second,  autlreflectlon  coatings  or  grooves  can  be  cut  into  the  silicon  section 
to  obtain  better  matching  from  metal  waveguide  to  silicon.  In  any  case,  we 
have  demonstrated  experimentally  that  there  is  very  little,  if  any,  elec¬ 
tronic  radiation  emanating  from  the  walls  of  the  exposed  silicon  guide.  This 
is  a  verification  of  the  calculations  which  indicate  only  a  very  small  frac¬ 
tion  of  power  is  located  outside  the  semiconductor  guide  vails,  and  that 
indeed  most  of  the  power  is  being  transmitted  inside  the  semiconductor. 

CONCLUSIONS 

An  analysis  was  made  of  waves  in  a  dielectric  slab.  In  considering 
this  structure,  it  was  expected  that  some  insight  would  be  obtained  relating 
to  the  fraction  of  power  propagating  inside  a  dielectric  waveguide. 

The  properties  of  the  slab  were  chosen  to  match  possible  experiments 
with  high  resistivity  silicon  at  frequencies  associated  with  Ku  band. 

That  Is,  In  the  calculations  chosen  it  was  assumed  that  silicon  has  a  di¬ 
electric  constant  of  12  and  can  be  obtained  lossless  (O'  ■  0).  The  height 
of  the  dielectric  slab  was  taken  at  the  same  dimensions  as  the  inner  height 
of  metallic  waveguide  for  Ku  band  (Figure  3  and  4),  and  a  TM  mode  was  cal¬ 
culated. 

For  the  TM  mode  in  the  dielectric  slab  with  materials  and  dimensions 
given  it  is  seen  that  three  possible  modes  can  be  naintained. 

Mode  1,  H°  =  cos  k,  y  ,  k  =  .263x10^ 

a  11 

Mode  2,  H°  =  sin  k  y  ,  k  =  .78x10^ 

*  1  1 

Mode  3,  IL  =  cos  k  y  or  k  =  1.09xlo3 

x  1  1 

If  we  address  ourselves  to  the  question  as  to  where  most  ex.  the 
energy  is  located,  we  find  that  In  all  three  cases,  most  of  the  power  flows 
In  the  dielectric.  This  is  a  significant  point,  since  if  we  want  to  con¬ 
struct  active  devices  in  the  semiconductor  dielectric,  it  is  necessary  for 
these  devices  to  be  in  a  location  where  they  can  interact  with  the  energy 
being  propagated. 
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In  addition,  we  find  that  for  the  lower  modes,  the  energy  tends  to  he 
confined  more  tightly  into  the  dielectric.  This  tendency  appears  as  indica¬ 
ted  in  the  following  figures.  (Figure  11). 

The  situation  now  appears  hopeful  for  guiding  the  electromagnetic  wave 

rrrriTM  ^TT.I  2Cr«3UCtCrt 

In  the  next  phase  of  work  we  must  analyze  the  mode  characteristics  for 
a  rectangular  waveguideB.  We  shall  need  this  information  in  order  to  know 
best  how  to  launch  this  wave  from  metal  waveguide  to  semiconductor  dielectric 
and  vice  versa  with  a  minimum  of  losses  or  reflections.  In  addition,  the 
coupling  from  the  semiconductor  guide  to  the  devices  will  be  of  prime  concern 
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GLOSSARY 


'Si.  “ 


■  frequency  =  15- ^  OH2 

=  wavelength  In  free  space  =  I.9U  cm 
=  wavelength  in  silicon  ---  .55  cm 


tT 


2a  =  cutoff  wavelength  =  3.16  cm 

free  space  Impedance  «  37^  ohms 

Impedance  In  silicon  =  107  ohms 

conductivity  =  5*^  x  10"2  (Am)  ^ 

attenuation  of  silicon  in  guide  =  2.89  neper/m 

p  ■  impedance  in  air  filled  waveguide  = 

*-oi 


576 

^  =  ^6  °hm 

Impedance  in  silicon  filled  guide 


4k  - 


v  / _ 

iT~Tm  1  r  -  y  v  1 3 T*  c  107  ohms 

rr 

propagation  constant  in  silicon  filled  guide 

ft OQ. 

wavelength  in  silicon  in  guide 


wavelength  in  air  filled  guide 


Block  length  =  10  cm  ~  18  wavelengths 

7?5B- 
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FIG,  5 

POSSIBLE  DISTRIBUTION  OF  FIELD  IN  DIELECTRIC 


I 


FIG. 7 

PENETRATION  OF  FIELD  OUTSIDE  SURFACE 
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1.603,-4.097 


1.6034,-4.097 


FIG .  8 

GRAPHICAL  SOLUTIONS 


FIG  9 

COSINE  MODE 


FIG.  10 

SECOND  COSINE  MODE 


TEST 

SECTION 


SCOPE 


ATTENUATOR  PHASE  DETECTOR 
SHIFTED 


FIG.  12  EXPERIMENTAL  BRIDGE  CIRCUIT 


CONDITION 

© 


TEST  SECTION 


4  1/2 


Y//////Z//////A  **- - SILICON  *2000  OHM-CM 


T 


T 


4  7/32" 


W///////AW, d 


1 

i 

I 

I 

f- — 2  1/2" 


i 

i 

I 

j 
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i 

i 

l 


FIG.  15  EXPERIMENTAL  ARRANGEMENTS 
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MODEL  FOR  ANALYSIS  OF  MULTIPLE  REFLECTIONS  IN  SILICON 


APPKtmy  a 


Calculation  of  Attenuation  Due  +.n  <?■(  n  +v«  i?™™4 _ 


The  attenuation  due  Lu  the  conductivity  and  the  multiple  reflections 
within  the  semiconductor  slat  used  in  these  experiments  can  be  calculated 
from  the  ratio  of  the  transmitted  electric  field  to  the  incident  electric 
field.  This  ratio  is  given  by 

—  =  Ytf{'c; ^  -*  -2^.  5/h4/7*/) ,  (i) 

V  Z*, b  / 

where  the  transmission  coefficient , 

>fr  -  J-&*  (2) 

Z-a,b  Zo, 

and  2„  ,  the  impedance  of  waveguide  in  air  is 


37  A 


rti' 


the  lmPedance  at  the  first  interface  of  the  dielectrics  is 


?  .  —  7Z  f 

— - — - -p 

n  °  3.  'fduhfep. 


t)' 


as  in  Figure  (lhA )  where  7*  V  —  rs 

and  C  =  ^  v4a-‘ y "  " C  ' 

The  impedance  in  the  semiconductor  is  approximately 


_  3ZL- 


rr 


(r  -  /£. 


Substituting  into  (3).  P*  _  7  /) 
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and  the  identities, 


tanh  (x  +  y)  =  '^‘Qhhx  V-  V 

I  n  h  X  fa.hk  y 


■'!  Oy  ■ 


equation  (3)  becomes 


-  2 ■ 

A~  j  7c5  *  4  °Ac<i>  JA / 


When/Ct/"  hTT. 


Substituting  into  (1)  -f^j_  -  °^e£. 

and  the  identities, 


,^L  0,9 


aJ 


Cosh  (tty)  ■-  Cosh  /  Co%k y  -+-  5/nhx  s/hhy 

5'*h(x+y)  -  5/nh  /  Ccsh y  +•  <A°5*X  S>*h  y 


equation  (1)  becomes 


-f“  -  •*„  ( _  ^2  S.nh  f\ 

1  / 

At  15.40  GHz  (  A  air  =  JL.94  cm)  the  10  cm  long  silicon  slab  is  a  half 
transformer  making  /%$  -r  =477  .  The  /\„  cutoff  vavelength  for  the 


metal  guide  is  3*16  cm,  therefore 


-  1/ - T^T - TX  =  oims 

(/-(&) 


Z  -  i  7i 


n  / 0  7  ohm  5 


^ 


^  £Z"  -  C'&Z'i  htysrs  yc 


r  tP-Cto  leiijtk  , 


-  Xvx/f 


hi  ffijtff 
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Inserting  these  values  into  equations  (3)  and  (2)  yields 

/  U  71.  -4-  !  ft  7  4-  .  O  Jt  \ 

/  _  /W  -  -  -  n  oU, 

*  \/°7  4-76  ft*  k  '??clJ 

Y  -  *(**>1-  -  ~  P.  6  y  . 

r  M-h?  76 


Equation  (l)  then  becomes, 

to  —  __  /£7  S/tilt.tff 


-  a  y# 


and  the  ratio  of  pover  transmitted  is 

l-~-i  -  H\  .33/  .r  </.  SdB. 


The  measured  value  of  attenuation  in  silicon  filled  waveguide  was 
approximately  4  dB.  Considering  the  range  of  resistivities  measured  over 

/,o-q  l0T^  Sl^c°n  ®labj  ’'1°5°  "  3300  ohm-cm)  the  average  value  used, 
(1830  ohm-cm)  could  be  in  error  by  approximately  20%.  An  agreement  of  cal- 
culated  vs  measured  attenuation  in  this  case  is  considered  valid 
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APPENDIX  B 


r  :  '-Q*  •  cO  Ocuati oru:  jC  ar.d  2-1 


Fi.  r.sOSc'.vA  -34/21/71  4;22  PM. 


■  ■  u  ;  AR  j  } 

0,;  „  ">  f  j  9  Y  " 

- r'j t  L3A  START2 
5  'J  it  3  CARD 
633  S i A  COT 
730  i  CARD 

*00  xSA  FI  iNID  ALLROO TS 
<500  F  CARD 

(000  PEAL  PROCEDURE  FCT(X);  REAL  X; 

1100  FCT: =X*2\ (l+COT  (X)*2/144)-l 9.36} 

1200  REAL  ARRAY  P,  Q,  CK{0:20]; 

1300  INTEGER  K; 

U00  FORMAT  FMTKX4,  "P",  XS,  "Q",  X8,  "CK")j 
1533  FINDALLROOTSC  FCT,  -4.4,  4.4,  100,  9,  p)  ; 
1^02  FOR  K:  =  i  STEP  1  UNTIL  P[ 0  ]  +  . 1  DO 
i70o  0[  K 1 ; -  —  P  C  K 1\C0T(P[K] )/12; 

!*0J  CKCKl!=PtK]\PIK]+QtK]\flIK 1-19.3 6; 

1900  WRITE (TYPE,  FMT 1 ) } 

2000  FOR  Xtsl  STEP  1  UNTIL  Pt 0  ]+. 1  DO 
2100  WRITE (TYPER,  PtKJ,  Q f  K  ] .  CK[K1): 

2200  $$A  FI N I 
2300  END. 


END  QUIKLST  .8  SEC. 


RUNNI NO 

TIME  ON  1603  21  APR  1971 


CK 


P  G 

-  1  .3 9*3 93+00-  1  . 1  90293-0  1 
*3  .23  .  :--3  3+00-2.986  1  73*00 
*3.06,063*00-  3.16068@*00 
0.061063+03'  3.160633+00 
3.23!  533+001-2 . 586  1  7(34-00 
. 3  0  *  39  3*0 0-  1  . 1  90 293-  0  1 

I  Ml-  1603 


0.000003+00 
0.000003+00 
2.000003+00 
3.020003+00 
0 . 002003+  00 
0 .00000  '3+  0  0 


CSEGoA  2.0  3 c.C  . 
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FILEjCSEQIIA  -04/22/71  10*35  AM. 


100 
200 
300 
400 
500 
600 
700 
800 
900 
1000 
1  100 
1200 
1300 
1400 
1500 
1600 
1  700 
1800 
1900 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
2800 
2900 
3000 


BEGIN 
S*A  STAPT1 
HENRY" 

S$A  START2 
$  CARD 
$$A  TAN 
$  CARD 

SSA  FINDALLROOTS 
£  CARD 

REAL  PROCEDURE  FCT<X)j  REAL 
FCTi =X*2\(1+TANCX>*2/144 
REAL  ARRAY  P  0,  CKt0*201* 
REAL  PI.  01; 

INTEGER  K{ 

FORMAT  FMT1(X4,  "p"  X8  "  (?' 
FI NDALLROOTSCFCT  -4.4  ’4.4 
FOR  Kt=l  STEP  1  UNTIL  P[0]+! 
CCK]x=P(Kl\TAN(PtK])/12i 
CKm,:P[K]\p[KHOfK]\On 
WRITE(TYPE,  FMT1 ) J 
FOR  Ks=1  STEP  1  UNTIL  Pr0)+. 

WRITECTYPER,  Pf k ] .  Q[ K3 
WRITE(TYPER) •  '  ’ 

FOR  N,  =.4. 4  STEP  .1  UNTIL 
BEGIN 

01l=Pl\TANCPl)/12x 
e»™ITE<TYPE*'  P1’  on? 

$$A  FINI 
END. 


END  OUIKLST  1.3  SEC. 
COMPILING. 


END  COMPILE  9.7  SEC. 
RUNNING 


TIME  ON  1036 


22  APR  1971 


X? 

>-19.36} 


,  X8,  CX  > ; 

100,  9,  P); 
1  DO 

1-19.36; 

1  DO 

CKtKD? 

4.41  DO 
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P  f  CK 

-4. 315529+00  8.580769-01  0.0  0  0  0  0  94-00 
-1  .603399+00-4  .097459+00  0.000009+00 
-1.539689+00  4.121829+00  0.0000004-00 
i  . 5396t<ii'-*-00  a  .  i2i  k .0010 0 kk*4- 100 

1  .603399+00-4 .097459+00  0.000009+00 
4  .31  552€4.00  8.580769-01  0.000009+00 


-4  .4000094-00  1.13532(24-00 
-4.300009+00  8.19095®- 01 
-4  .200009+00  6.222239-01 
-4.1  0000®4-00  4  .86372(3-0 1 
-4  ,00000®4-00  3 . 85940®-0  1 
-3  .900009+00  3. 07913®- 01 
-3  ,80000®4-00  2  .44959@-0 1 
-3.70000®4-00  1  .926269-01 
-3.600009+00  1  «48040®-0 1 
-3.50000®4-00  1  .09254@-0  1 
-3.40000S+00  7.48898®-02 
-3  .300009+00  4 .3930 1@-02 
-3  ,20000®4-00  1  .559300-02 
-3,1 03009+00-  1  .0751  0®-02 
-3  . 0000004-00-3. 563669-02 
-2.900009+00-5.954809-02 
-2  .80000®4-00-8  .29570@-02 
-2 , 70000®4-00- 1  ,06364®-0l 
-2.600009+00-1  .303460-01 
-2.500009+00-1 .55630®-01 
-2.400009+00-1 .832039-01 
-2 ,30000®4-00-2 .14  51  6@~01 
-2  .20000SH-00-2. 51  868®-01 
-2 .1000 09+00-2. 99223®- 01 
-2  .00000®4-00-3 .64 1  73®-01 
-1  ,90000®4-00-4.63457®-01 
-1 . 800009+00-6. 42939®-01 
-1  , 70000®4-00- 1  »09035®+00 
-1 .60000©+00-4.56434®+00 
-1.500009+00  1.76268®+00 
-1.400009+00  6.76420®-01 
-  1  .300009+00  3 . 902289-01 
-1.200009+00  2r57215@-01 
-1.1000094-00  1,801  03®-0 1 
-1.000009+00  1.297849-01 
-9.000009-0 1  9.451199-02 
-8.000009-Cl  6.864269-02 
-7.000009-01  4.913359-02 
-6.000009-01  3.420689-02 
-5.000009-01  2.276269-02 
-4.000009-01  1.409319-02 
-3.000009-01  7.733419-03 
-2.000009-01  3.378509-03 
-1.000009-01  8.361229-04 
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2. 009989-10  3.366699-21 
1.000009-01  8.361229-04 
2.000009-01  3.378509-03 
3.000009-01  7.733419-03 
4.000009-01  1.409319-02 
5.000009-01  2.276269-02 
6,000009-01  3.420689-02 
7.000009-01  4.913359-02 
8.000009-01  6.864269-02 
9.000009-01  9.451199-02 
1.000009+00  1.297849-01 
1 . 1 00009+00  1.801039-01 
1.200009+00  2.572159-01 
1.300009+00  3.902289-01 
1.400009+00  6.764209-01 
1.500009+00  1 .762689+00 
1 .600009+00-4.564349+00 
1 .700009+00-1 .090359+00 
1 .800009+00-6.429399-01 
1 .900009+00-4.634579-01 
2 .000009+00-3 . 64 1 739-0 1 
2 . 1 00009+00-2.992239-0 1 
2  .200009+00-2.518689-01 
2 .300009+00-2 .1451 69-01 
2.400009+00-1 .832039-01 
2.500009+00-1 .556309-01 
2.600009+00-1 .303469-01 
2.700009+00-1 .063649-01 
2 .800009+00-8 .295709-02 
2 .900009+00-5.954809-02 
3.000009+00-3.563669-02 
3.100009+00-1.075109-02 
3.200009+00  1.559309-02 
3.300009+00  4.393019-02 
3.400009+00  7.488989-02 
3.500009+00  1.092549-01 
3.600009+00  1.480409-01 
3.700009+00  1.926269-01 
3.800009+00  2.449599-01 
3.900009+00  3.079139-01 
4.000009+00  3.859409-01 
4.100009+00  4.863729-01 
4.200009+00  6.222239-01 
4.300009+00  8.190959-01 
4.400009+00  1.135329+00 

TIME  OFF  1037 


END  CSEQUA  3.7  SEC. 
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